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Abstract. We discuss a novel dark matter signature relevant for directional detection of 
Weakly Interacting Massive Particles (WIMPs). For heavy enough WIMPs and low enough 
recoil energies, the maximum of the recoil rate is not in the direction of the average WIMP 
arrival direction but in a ring around it at an angular radius that increases with the WIMP 
mass and can approach 90° at very low energies. The ring is easier to detect for smaller 
WIMP velocity dispersion and larger average WIMP velocities relative to the detector. In 
principle the ring could be used as an additional indication of the WIMP mass range. 
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1 Introduction 

Direct dark matter experiments search for energy deposited in low-background detectors by 
the scattering of Weakly Interacting Massive Particles (WIMPs) present in the dark halo of 
our galaxy. If a signal is observed, it is necessary to have unmistakable signatures that it is 
due to dark matter. Two such signatures are the annual modulation of the signal due to the 
motion of the Earth around the Sun, extensively studied since [1] , and the daily modulation 
due to the spinning of the Earth around its axis. 

With the advent of directional detectors [2] using CS2, CF4 or 3 He, such as DRIFT [3], 
DMTPC [4], NEWAGE [5] and MIMAC [6], the observation of other dark matter signatures 
becomes possible. Because dark matter WIMPs arrive to us from a preferential direction, one 
expects an anisotropy in the recoil event rate [7] . The detector is moving with a velocity Vi a b 
with respect to the Galaxy. If the particles in the dark halo of our galaxy are on average at 
rest with respect to the Galaxy, the average velocity of the dark matter particles with respect 
to the detector is — Vi a b. This is the case in the standard halo model. Here we consider a 
galactic dark matter halo model whose velocity distribution, as in the standard halo model, 
is locally approximated by an isotropic Maxwell-Boltzmann (1MB) distribution, but with a 
velocity dispersion independent of the average velocity, to better represent what is known 
about the actual local parameters of the distribution. We use the 1MB in most of this paper. 

The easiest dark matter signature to detect in directional detectors is the average recoil 
direction, expected to be in the direction of — V"i a b [7-10]. Ref. [10] finds that for an S 
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Figure 1. Figure describing the existence of a ring in recoil q space due to a WIMP, or WIMP 
stream, with velocity v. \i is the reduced WIMP-nucleus mass. 



detector with an energy threshold of 20 keV, 9 events would be enough to reject isotropy 
at 95% CL, while between 27 and 32 events, depending on the velocity distribution, would 
be enough to confirm the direction of solar motion as the median inverse recoil direction at 
95% CL. For a 100 GeV/c 2 WIMP, with spin-independent interactions and WIMP-proton 
cross section a p = 10~ 44 cm 2 , 30 events in CS2 corresponds to an exposure of 150 kg-yr [9]. 
Ref. [10] assumes the recoil directions, including their senses, can be reconstructed perfectly 
in 3d and the background is zero. Assuming a less optimistic model for a CF4 detector with 
a 50% background contamination, Ref. [11] finds that a signal pointing within 20° of the 
opposite direction to the Solar motion can be confirmed with 25 events between 5 keV and 
50 keV. 

In this paper we discuss an additional dark matter signature relevant for directional 
detection as a secondary feature with respect to the mean recoil direction: a ring of maximum 
recoil rate around the direction of — Vi a b- For heavy enough WIMPs and low enough energies, 
the maximum of the recoil rate is not in the direction of the average WIMP arrival direction 
but in a ring around it at an angle which can be up to 90°. 

The origin of the ring can be understood considering just one incoming WIMP with 
velocity v. Energy and momentum conservation in the collision of this WIMP with a nucleus 
of mass M imply that the magnitude of the recoil nuclear momentum is q = 2jiv cos 9, where 
v = |v|, \i is the reduced WIMP-nucleus mass and 9 is the scattering angle, i.e. the angle 
between v and the recoil momentum q [12]. As shown in Fig. 1, in momentum space these 
values of q lie on the surface of a sphere of radius fj,v, centered at /xv and passing through 
the origin. For a certain recoil energy Er = q 2 /2M, q must be, in momentum space, on 
a sphere with radius q = \J2MEr and centered at the origin q = 0. The only possible q 
values at fixed magnitude q are therefore on the ring on which the two spheres mentioned 
intersect. The angular radius of the ring 7 fulfills the condition cos 7 = \/2MEr/2[iv. In 
the 1MB, most WIMPs arrive with velocity — Vi a b, so the maximum number of recoils occur 
when v = Vj a b (see Eq. 2.5). As V| a b increases, the angle 7 increases, and the flux in the 
direction of — Vi a b decreases with respect to the maximum rate. This makes the ring easier 
to observe. 

Our ring-like signature would require more events than those mentioned above to be 
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detected. In Section 3 we find with a simple statistical test that for a CS2 detector with 
3d track reconstruction and a 100 GeV/c 2 WIMP, for recoil energies between 5 keV and 40 
keV and no background, 41 events total over the whole sky would be required to observe 
the ring in one of the most favorable cases (not taking into account the energy and angular 
experimental resolutions). This corresponds to an exposure of 250 kg-yr for a p = 10 -44 cm 2 . 
For recoil energies between 20 keV and 40 keV instead, 55 events in total over the whole sky 
would be required, corresponding to 700 kg-yr for a p = 10 -44 cm 2 . These are lower bounds 
on the number of events needed to detect the ring. 

We noticed the existence of this ring before (see e.g. Fig. 4 of Ref. [13] and the explana- 
tion of this figure) but to our knowledge its potential relevance for directional detection has 
not been pointed out. 



2 Differential Recoil Spectrum 



It is easy to see that at low recoil energies for heavy enough WIMPs the maximum of the 
directional recoil rate is not in the direction of — Vi a b, by writing the recoil spectrum in terms 
of the Radon transform. Let us start with the directional differential recoil spectrum as a 
function of the recoil momentum q [14] 
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Here p is the dark matter density in the solar neighborhood and m is the WIMP mass. 
The sum is over the nuclear species i in the target, and C% and ^ are the mass fraction 
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and the reduced WIMP-nucleus mass for nuclide i, respectively. dVt q 
an infinitesimal solid angle around the recoil direction q = q/g, q = 
of the recoil momentum, q/2fi = v q is the minimum velocity a WIMP must have to impart 
a recoil momentum q to the nucleus of mass M, or equivalently to deposit a recoil energy 
Er = q 2 /2M, /i = mM/(m + M), and /i a b is the 3-dimensional Radon transform of the 
WIMP velocity distribution. (Ji{q) is the WIMP-nucleus scattering cross section which can 
be split into spin-independent (SI) and spin-dependent (SD) parts, c?i(q) = crf l (q) + af T) (q). 

In the usual 1MB model, WIMPs are on average at rest with respect to the Galaxy, and 
have a Maxwellian velocity distribution with dispersion a v , truncated at the escape speed 
v esc (with respect to the Galaxy). Normalized to 1, the WIMP velocity distribution in the 
laboratory rest-frame is given by [12] 
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The laboratory is moving with velocity Vi a b with respect to the Galaxy (thus — Vi a b is the 
average velocity of the WIMPs with respect to the detector). Following Ref. [15], which gives 
100 km/s and 130 km/s as extreme estimates for the ID velocity dispersion a v /\^3, we take 
a v either 173 km/s or 225 km/s. A recent study by the RAVE survey using high velocity 
stars finds an escape speed in the range 498 km/s < v csc < 608 km/s at 90% confidence, with 
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Date 


VbaiRot (km/s) 


|V G alRot + V So i ar | (klll/s) 


Viab (km/s) 


December 1 


180 


193 


179.8 


May 30 


180 


193 


208.8 


December 5 


312 


324 


310.5 


June 2 


312 


324 


340.2 



Table 1. Maximum and minimum values of Vjab in km/s, for our two choices of the Galactic rotation 
speed VbaiRot ■ 



a median-likelihood value of 544 km/s [16]. We use the median RAVE value, v csc — 544 km/s 
for most of our results, but for a few figures we keep the older usual value t> esc = 650 km/s. 
For the Galactic rotation speed VcaiRot at the position of the Sun we take 180 km/s and 312 
km/s as low and high estimates (see Appendix A). 

Due to the ellipticity of the Earth's orbit the times of maximum and minimum Vi a b are 
not exactly half a year apart. These times depend on VcaiRot- For VcaiRot = 180 km/s, VJ a b is 
maximum on May 30 and minimum on December 1. For VcaiRot = 312 km/s, the maximum 
and minimum happen on June 2 and December 5, respectively. Table 1 gives the maximum 
and minimum values of |Vi a b| = Mab (in km/s) for our two choices of VfcalRotj including the 
contribution of the Solar motion Vg i ar . 

The Radon transform in the laboratory frame for the truncated Maxwellian WIMP 
velocity distribution in Eq. 2.2 is [12] 
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if (q/2/j,) + q • V"i a b < v esc , and zero otherwise. 

The recoil momentum q is measured in a reference frame fixed to the detector. The 
detector frame is at some orientation in the laboratory frame, which we define as fixed to 
the Earth with axes pointing to the North, the West and the Zenith. The transformation 
equations between the detector frame and the laboratory frame can be conveniently written 
in terms of direction cosines measurable in any experiment. This is formulated in Appendix 
A, where we also give the transformations from the laboratory frame in any location on 
Earth to the Galactic reference frame. The transformations in Appendix A take into account 
Earth's rotation around its axis, which is usually neglected. 

From Eq. 2.4 we can immediately see that if q/2^ < |Vi a b|, i.e. v q < Vj a b, the maximum 
of /lab (^ 5 q) occurs when — q • V[ a b = v q . Only if v q > Vi a b the maximum is in the direction 
of — Vi a b. Thus when the minimum WIMP speed v q required to cause a recoil momentum q 
is less than VJ a b, the maximum of /i a b occurs at an angle 7 between q and — Vi a b given by 
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The ring is clearly visible when /i a b is plotted on the sphere of recoil directions, as for 
example on panel (a) of Fig. 2. In this figure we use the HEALPix pixelization [17] to plot 
/lab on the sphere of recoil directions for sulfur recoils at two different times in the year and 
recoil energies, assuming VcaiRot = 280 km/s, a v = 173 km/s, v esc = 650 km/s, and m = 100 
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Figure 2. Angular recoil distribution /i a b(fq,q) at fixed v q for sulfur recoils of 100 GeV/c 2 WIMPs 
on (a) December 4 with Er = 5 keV, thus v q = 113 km/s, and (b) June 2 with En = 40 keV, thus 
v q = 319 km/s. The distribution is plotted on the sphere of recoil directions q using the HEALPix 
pixelization. The 1MB Model is used with parameters VcaiRot = 280 km/s, a v = 173 km/s, and 
v csc = 650 km/s. The North and South Galactic Poles (NGP and SGP), the North and South 
celestial poles (N and S), and the Earth meridians and parallels are indicated. The arrows show the 
direction of the average WIMP velocity — Vi a b, directed into the South Galactic Hemisphere (SGH) in 
December and into the North Galactic Hemisphere (NGH) in June. The color scale/grayscale shown 
in the vertical bars correspond to equal steps between the minimum and maximum values in units of 
(km/s) . Notice the ring of maximum /i a b values in (a) at 7 = 66° of — Vi a b- 

GeV/c 2 . Fig. 2. a shows /i a b on December 4 (when V[ a b is minimum) for En = 5 keV. We can 
see a ring of maximum /i a b, i.e. of maximum recoil directional rate around — Vi a b- — V] a b 
points slightly towards the South Galactic Hemisphere. Fig. 2.b shows /i a b on June 2 (when 
Vi a b is maximum) and for Er = 40 keV. In this case the maximum of /i a b and thus of the 
recoil rate, is in the direction of — Vi a b- — Vi a b points slightly towards the North Galactic 
Hemisphere at this time. The color scale/grayscale plotted on the spheres in Fig. 2 indicate 
different values between 0.00018 (km/s) 1 and 0.0023 (km/s) -1 for Fig. 2. a and between 
1.3 x 10~ 6 (km/s) -1 and 0.0023 (km/s) -1 for Fig. 2.b as shown in the horizontal bar. The 
minimum WIMP speed is v q = 113 km/s and Vi a b = 278.7 km/s in Fig. 2. a and v q = 319 
km/s and Vj a b = 308.3 km/s in Fig. 2.b. 

From Eq. 2.5, the ring of maximum /i a b exists if Er < 2{j?V^ h /M. For example, for 
WIMPs of mass m = 100 GeV/c 2 and V[ a b = 278.7 km/s, as in Fig. 2. a, a ring exists in S 
recoils if Er < 30.6 keV. Since the reduced mass approaches its maximal value [i = M as 
m s> M, there can be no ring if En > 2MV^ h for any value of m and the maximum possible 
ring radius (and thus the minimum possible cos 7) is given by 



For Er = 2MVj 2 b , 7 max becomes zero. Fig. 3 shows plots of 7 max as a function of Er when 
Viab is maximum for different target elements relevant to directional detectors. Here, either 
VbaiRot = 180 km/s or 312 km/s. The allowed values of 7 are below the lines shown. 
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Figure 3. Maximum ring radius 7 max , measured in degrees from the direction of — Vi a b, as a function 
of En, for three target elements S, F, and C. The labels indicate the target element, the subscripts 
showing the value of VoaiRot in km/s. The sulfur solid black line Siso coincides with the carbon 
dashed orange line C312. The possible values of the ring radius 7 lie below the respective curves, 
and to the left of the vertical lines (labeled E$, Ep, and Ec for S, F, and C, respectively), where 
/center//ring < 0.9 (see Eq. 2.11). The ring would be very difficult to observe with a smaller contrast, 
namely for larger /c e „t C r//rin g - 



The ring- like recoil feature we discuss here depends only on /i a b- In particular, it is 
independent of the type of WIMP interaction. However, to give concrete examples of how the 
ring would appear in an actual recoil rate, we show the ring for WIMPs with spin-independent 
interaction and equal couplings to protons and neutrons in a CS2 detector. For these WIMPs, 
the cross section in Eq. 2.1 is af l (q) = pf Afa p Si(q) / where p p = mm p /(m + m p ) is the 
WIMP-proton reduced mass, a p is the WIMP-proton cross section, Ai is the mass number 
of the nuclear species i, and Si(q) is the nuclear form factor, for which we use the Helm [18] 
nuclear form factor normalized to 1. The SI directional recoil rate is therefore [14] 



dR 
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Here po.3 is the dark matter density in units of 0.3 GeV/c 2 /cm 3 , 044 is the WIMP-proton 
cross section in units of 10~ 44 cm 2 , p p and m are in GeV/c 2 , and /i a b is in (km/s) -1 . The 
standard value of the local dark matter density is p = 0.3 GeV/c 2 /cm 3 , and this is what we 
use here. However, one should keep in mind the large uncertainties in this parameter. Recent 
astronomical constraints are consistent with 0.2 GeV/c 2 /cm 3 < p < 0.4 GeV/c 2 /cm 3 for a 
spherical dark matter halo profile, and up to 20% larger for non-spherical haloes [19]. Using 
a halo model independent method, Ref. [20] finds p = 0.43 ± 0.11 ± 0.10 GeV/c 2 /cm 3 , with 
uncertainties from two different sources. Ref. [21] finds the density for a specific simulated 
galaxy resembling the Milky Way is 21% larger than the mean value of p = 0.39 GeV/c 2 /cm 3 
obtained in a previous study [22] in which spherical symmetry was assumed. 

Fig. 2. a shows the ring at a particular time. However, in direct dark matter searches 
the event rate will need to be measured during several years to obtain the required exposure 
to observe the ring. Fig. 4 shows that the recoil rate averaged over a year is very similar to 
the distribution computed at a particular time. The figure shows the plot of the directional 
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Figure 4. Directional differential recoil rate as a function of the polar angle 9y measured from — Vi a b 
for VcaiRot = 312 km/s, a v = 173 km/s, w csc = 544 km/s, m = 100 GeV/c 2 , E R = 5 kcV in CS 2 . 
The solid black curve is the average rate over the year 2011, and the dashed green curve is the rate 
on June 2, when Vj a b is maximum. 



rate in CS2 averaged over one year (solid black curve) and on June 2 (dashed green curve) 
as a function of the opening angle 6y between the recoil direction and the average — Vi a b 
over the year, at one fixed azimuthal angle (j>y. In the 1MB, for a given 9y the rate is the 
same at different azimuthal angles around — Vi ab . Here VcaiRot = 312 km/s, a v = 173 km/s, 
v C sc = 544 km/s, m = 100 GeV/c 2 , and Er = 5 keV. If a ring exists, the rate reaches its 
maximum at 0y = 7. The two curves shown in Fig. 4 are very similar, thus the effect of 
integrating the rate over time is not significant. For simplicity, in this paper we compute the 
rate at specific dates only (we choose those at which Vj a b is maximum). 

The plots of /iab or directional rate are better viewed not in a HEALPix sphere, but in 
a planar projection. In Fig. 5 and 6 we use the Mollweide equal-area projection maps of the 
celestial sphere in Galactic coordinates. The relationship between the (x,y) coordinates on 
the Mollweide map and the Galactic longitude and latitude (l,b) is given by [23] 

-ttx , /20 + sin(2(9)\ , . 

I = — = , b = arcsin , (2.8 

2V2cos0 V 7T /' v ' 

where 

6 = arcsin (y/V2). (2.9) 

We first show, in Fig. 5. a, a Mollweide projection of the number fraction of WIMPs 
-Fwimp(v, v q ) crossing the detector per unit solid angle with speed larger than function 
of the WIMP velocity direction v (these WIMPs produce a recoil momentum of magnitude 
q or higher when scattering off a nucleus of mass M) [13] 

f»max(v) 

FmMp(v,v q )= /wiMp(v)-y 2 dv. (2.10) 

Jv„ 
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Figure 5. Mollweide equal-area projection maps of the celestial sphere in Galactic coordinates 
showing (a) the number fraction i r wiMp(v,w 9 ) of m = 100 GeV/c 2 WIMPs crossing the Earth per 
unit solid angle as a function of the WIMP velocity direction v. For this figure we take a minimum 
speed v q = 113 km/s, as necessary to produce En = 5 keV sulfur recoils, (b) The directional 
differential recoil rate in CS2 at En = 5 keV for m = 100 GeV/c 2 . In both panels we assume the 
1MB with v esc = 544 km/s, a v = 173 km/s and VcaiRot = 312 km/s on June 2. Notice the direction 
of — Vi a b marked with a cross. The color scale/grayscale shown in the vertical bars corresponds to 
equal steps between the minimum and maximum values in 5. a in units of sr _1 , and in 5.b in units of 
10~ 6 x (po.soWkg-day-keV-sr). Eq. 2.11 gives /center //ring = 0.42 in the right panel. 



The upper limit of this integral is w m ax(v) = -v • V lab + y (v • V lab ) 2 - Vf ab + v% sc and 
the analytic expression of i*wiMp(v, u<j) is given in Eq. 13 of Ref. [13]. The maximum of 
^WIMp(v, v q ) happens when v • V^b = — Vjab> i-e. in the direction of the average WIMP 
velocity — Vi a b- Most WIMPs move in the direction opposite to the laboratory motion, 
marked by a cross in the figures. 

In Fig. 5.b we show a Mollweide map of the directional differential recoil rate in CS2, 
Eq. 2.7, produced by the WIMPs in Fig. 5. a in which the ring of maximum rate around the 
— Vi a b direction is clearly visible. In Fig. 5, we used m = 100 GeV/c 2 , and the 1MB with 
v C sc = 544 km/s, VcaiRot = 312 km/s and a v = 173 km/s on June 2. In Fig. 5.b the recoil 
energy is Er = 5 keV. 

It is easier to see the ring when the contrast between the rate at the center of the ring 
(in the direction of — Vi a b) and the ring is larger. In terms of /i a b, the ratio of the value 
/center at the center of the ring to the value / r i ng at the ring is approximately, for the 1MB 
neglecting the escape speed, 
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/center//ring = 0.42 in Fig. 5.b (see the 5 keV profile in Fig. 7.a). The smaller the ratio 
/center //ring 1 the easier it is to detect the ring. Thus the best prospects to observe the 
ring are at low recoil energies and for heavier WIMPs (so v q is small), large Vi a b and small 
a v . In Figs. 5.b and 6 we vary V^b and a v giving the four combinations of maximum and 
minimum values for both. Fig. 6 shows plots of the CS2 directional rate for Er = 5 keV, 
m = 100 GeV/c 2 , and three combinations of VJ a b and a v different from those in Fig. 5: (a) 
VcaiRot = 180 km/s, a v = 225 km/s on May 30; (b) VcaiRot = 180 km/s, a v = 173 km/s on 
May 30; and (c) VcaiRot = 312 km/s, a v = 225 km/s on June 2. The right panel of Fig. 5 
displays the fourth combination of VcaiRot and a v . The rate is dominated by scattering off 
S, by a factor of about 100. In Fig. 6. a, V^b = 208.8 km/s and a v = 225 km/s is the worst 
combination of low Vj a b and high a v . It is clearly seen from the figures that the ring is most 
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Figure 6. Directional differential recoil rate for different values of the Galactic rotation velocity 
VbaiRot and of the WIMP velocity dispersion a v : (a) VcaiRot = 180 km/s, a v = 225 km/s on May 30; 
(b) VcaiRot = 180 km/s, <r v = 173 km/s on May 30; (c) VbaiRot = 312 km/s, <r v = 225 km/s on June 
2. The right panel of Fig. 5 displays the fourth combination VcaiRot = 312 kms/ and a v — 173 km/s. 
All panels have En = 5 kcV in CS2, m = 100 GcV/c 2 , and v csc = 544 km/s (the same parameters as 
in the right panel of Fig. 5). The color scale/grayscale shown in the vertical bars correspond to equal 
steps between the minimum and maximum values given in units of 10~ 6 x (/9o.3C44/kg-day-keV-sr). 
In panels (a), (b), and (c), /center //ring = 0.91, 0.86, and 0.60, respectively. 



visible in Fig. 5.b, with the largest VJ a b and smallest a v combination of the four. This best 
combination is chosen for Fig. 7. 

The dependence of the ratio / center / J ring on energy can be extracted from Fig. 7. Fig. 7 
shows the CS2 directional differential rate as a function of 9y at different recoil energies for 
V| a b = 312 km/s and a v = 173 km/s on June 2. In Fig. 7. a, m = 100 GeV/c 2 , and the ring 
is present at all energies below 40 keV. Fig. 7.b shows the case of m = 600 GeV/c 2 at higher 
energies of E R = 40 keV, 50 keV and 60 keV. The ring is present at all energies below 50 
keV in this case. In Fig. 7.a, /center //ring = 0.42, 0.58, 0.80, and 0.99 for E R = 5, 10, 20, and 
40 keV, respectively. In Fig. 7.b, /center //ring = 0.89, 0.96, and 0.99 for E R = 40, 50, and 60 
keV, respectively. For the same energy, v q is smaller for a larger m. Thus as m increases, the 
ratio /center/Zring becomes smaller and the contrast between the center and the ring becomes 
larger. This is the reason that the ring is present for E R = 40 keV in Fig. 7.b but not in 
Fig. 7.a. 

For other materials the ring's dependence on energy is explored in Fig. 8, which shows 
the directional differential recoil rate for S, F, and C as a function of 9y at different recoil 
energies, 5, 10, 20, and 40 keV. In this figure again m = 100 GeV/c 2 , but we use the average 
values of a v and V[ a b, cr v = 199 km/s and VcaiRot = 246 km/s, which yields the maximum 
V| a b = 274.4 km/s on June 1, instead of the most favorable combination of values used in 
Fig. 7. Notice that since in CS2 the rate is dominated by scattering off S, Figs. 7. a and 8. a 
only differ because Vj a b is larger and a v is smaller, thus for energies below 40 keV, the ratio 
/center//ring is smaller and the contrast between the ring and the center is larger in Fig. 7. a. 

In Fig. 8. a, / CC ntcr//ring = 0.72, 0.85, and 0.97 for E R = 5, 10, 20 keV, respectively. In 
Fig. 8.b, /center/Zring = 0.78, 0.91, and ~ 1 for E R = 5, 10, and 20 keV, respectively. In both 
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Figure 7. Directional differential recoil rate as a function of the polar angle 9v measured from — Vi a b 
for CS 2 , with F Ga raot = 312 km/s, a v = 173 km/s on June 2, for (a) m = 100 GcV/c 2 at E R = 5, 10, 
20, and 40 keV (Fig. 5.b corresponds to the 5 kcV profile) and (b) m = 600 GcV/c 2 at E R = 40, 50, 
and 60 kcV. 




Figure 8. Same as Fig. 7. a but for (a) S, (b) F, and (c) C recoils in the 1MB with v csc = 544 km/s, 
cr.„ = 199 km/s and VbalRot = 246 km/s on June 1. Four recoil energies are shown in each plot: 5, 10, 
20, and 40 keV. 



Figs. 8.a and 8.b, v q > V iah for E R = 40 keV. In Fig. 8.c, /center //ring = 0.84 and 0.97 for 
En = 5 and 10 keV, respectively, and v q > Vi a b for Er = 20 and 40 keV. 

Fig. 8 illustrates the impact of using different targets on the existence of the ring in. 
When m 3> M, which is the case for all three panels in Fig. 8, v q ~ y Er/2M, and as the 
target mass M increases, v q decreases and the contrast between the ring and center becomes 
larger. In Fig. 8. a and Fig. 8.b, the ring is present at energies below 20 keV, whereas in 
Fig. 8.c, because C is lighter than S and F, the ring is present only at energies below 10 keV. 
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Figure 9. (a) Energy-integrated directional rate 1Z(Ei,E2,0v) (Eq. 2.12) in CS2 for different energy 
intervals as a function of 9y using the same parameters of Fig. 7. (b) Same rates re-scaled to be 
between and 1. 
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Figure 10. Same as Fig. 9. a but for 10 kcV energy intervals in (a) CS2 and (b) F. The black, 
green, orange (only in (a)) and red curves are for 10 keV < Er < 20 keV, 30 keV < Er < 40 keV, 
40 keV < E R < 50 keV and 50 keV < E R < 60 keV, respectively. 



Measuring the differential rate will require very large statistics. Thus we explore next 
the probability of observing the ring in the rate integrated over energy. 
Fig. 9. a shows the CS2 energy-integrated directional rate, 



£2 



dR 



n(E 1 ,E 2 ;6 v )= / JTP j- J 5- , dE R = - V' , , (2.12) 



El dEf{ d(j)y d COS 9y 



dR(E 1 ,E 2 ) 



d(j)y d COS Oy 



at fixed azimuthal angle (j)y integrated over different energy intervals as a function of Oy for 
m = 100 GeV/c 2 , Vi^ = 312 km/s and a v = 173 km/s on June 2. In Fig. 9.b, the integrated 
recoil rate has been re-scaled to be between and 1, what makes it easier to see that rings 
exist in the rate integrated from 5 keV up to 60 keV and lower energies. Fig. 10 shows 1Z 
for different 10 keV energy intervals instead, always as a function of Oy and for the same 
parameters of Fig. 9, again for CS2 (left panel), and also for F (right panel). Fig. 10. b is 
similar to Fig. 2 of Ref. [11] (in which no ring feature is seen), but for different values of Vj a b 
and a v . 

From Figs. 9 and 10 we can conclude that a ring can be observed in the integrated rate, 
but the energy threshold should be as low as possible. 
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3 Statistical significance of the ring 



Here we use a simple statistical test to estimate the number of events needed to detect the ring 
assuming no background (and not taking into account the energy and angular experimental 
resolutions). We define two angular regions centered at — Vi a b: the "ring" and the "center". 
The "ring" is defined as the region with 71 < By < 72 . The "center" is defined as the By < 71 
region. The angles 71 and 72, which determine the thickness of the ring, are in principle to 
be chosen in an optimal way. Here, for simplicity, we take them to be the By values for which 
the energy-integrated directional rate (at fixed azimuthal angle (f)y), 1Z(Ei,E 2 ; By) (Eq. 2.12) 
is halfway between its value at — Vj a b (i.e. By = 0) and its maximum value (at By = 7), 

niEu ft; 70 - n (El , g, ; 72 ) - ^ £ * °> + . (3.1) 

The solid angle subtended by the ring and the center regions so defined are Q r = 2-7r(cos7i — 
cos 72) and Q c = 2tt(1 — cos 71). 

To observe the ring at the 3a level, the difference between the number of events in the 
ring and the center regions, N r and N c , should be large enough. At the 3a level, then 

N r - N c > 3a ~ 3y/N, (3.2) 

where N = N r + N c is the total number of events for By < 72 over the duration of the exper- 
iment. The number of events needed depends on the energy interval. Assuming a particular 
interval (Ei,E 2 ), we approximate N r and N c using the middle values of 1Z (Eq. 2.12) in each 
region, 

N r (E x ,E 2 ) = ^ 7l ) +^ 72 ) o r MT = K r {E l ,E 2 ) fl r MT, (3.3) 

N C (E X ,E 2 ) = ^(°) + K-(Ti) Qc MT = K c (E h E 2 ) fl c MT. (3.4) 

MT is the detector exposure. To exemplify this procedure, Fig. 11 shows the ring thickness 
limited by 71 and 72, and the rates TZ T i Qg and 7£ cen t e r with dashed horizontal lines in the 
integrated directional recoil rate in the 5-40 keV interval reproduced from Fig. 9. a. 

From Eq. 3.2 to 3.4 we can eliminate the exposure MT (which depends on a p and p) 
and find N for the particular (Ei,E 2 ) energy range (for By < 72), 

Notice that the larger the difference in the solid angles fi r and £l c , the smaller is the needed 
-^min- Using the minimum number of events N m [ n (Ei, E 2 ) and Eqs. 3.3 and 3.4, we can find 
the minimum necessary exposure (MT) m [ n in units of kg-yr/po.3044) which allows us to find 
the minimum number of events in all directions iVt t (in any energy interval (E±,E 2 )), 

N tot (E[,E' 2 ) = (MT) min 2tt J Tl(E[, E' 2 ; By) d cos By = (MT) mill K tot (E[, E' 2 ). (3.6) 

Notice that iVtot is independent of a p p, as it can also be seen from an alternative expression 
for it, which avoids using the exposure, 

Nt ° tiE ^ E * ) - N {-Jl tot (E[,E>) + K^E') ) ■ (3 - ?) 
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Figure 11. Ring radius at By = 7 (central red vertical line) and thickness of the "ring" region 
(between the two lateral green vertical lines), which in this case extends from 71 = 26° to 72 = 62°, 
for the energy integrated directional rate in the 5-40 keV interval reproduced from Fig. 9. a. The 
dashed blue horizontal lines in the "central" and "ring" regions specify 7vL r i, lg , and lZ ccn t cr (Eqs. 3.3 
and 3.4), respectively. 



For CS 2 and m = 100 GeV/c 2 , V Uh = 312 km/s and a v = 173 km/s, in the 5-40 keV 
interval (see curve in Fig. 9), we find N m i n = 21 events for 9y < 62°, i.e. iV tot = 41 events 
over all directions (7£ ccn t r /Turing = 0.92, MT = 250 kg-yr/po.3044) are needed to detect the 
ring. In the 20-40 keV interval we would instead need iV tot = 55 events (T^centcr /Turing = 0.97, 
MT = 700 kg-yr/po.3C"44)- For this WIMP mass, the ring disappears when integrating over 
energies above 40 keV. With this last exposure, there would be 136 events above 20 keV (but 
the ring could not be seen if integrating above 40 keV). This is about 5 times more than 
needed to detect the mean recoil direction, according to Ref. [10]. 

With our simple method we can also estimate the number of events required to detect 
the forward-backward asymmetry Np-^, using an equation similar to Eq. 3.5. In this case, 
71 = 7r/2, and 72 = ir, thus O forward = ^backward = 27T. 7£ forward and ^backward are defined 
here as 

- ft(0)+ft(7T/2) 
'Morward(^l) &2) ~ „ ' 

^backward ( E\, Ey) = — ^ ^ ^ , (3.8) 

where the rates 1Z are also defined between E\ and E2 (Eq. 2.12). Thus 

N FB > U gforward + gbackward 'j (3 9) 

V ^-forward ^-backward / 

For CS2 and m = 100 GeV/c 2 , we find that iVpB = 17 events would be required 
above 5 keV to detect the forward-backward asymmetry ("/^backward /^-forward = 0.16, and the 
exposure needed would be 57 kg-yr/po.3044)- For Er > 20 keV instead, we find iYpB = 13 
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events (/^-backward /"^-forward = 0.09, MT = 56 kg-yr//9o.3C44) • This number, obtained with 
our simple method, compares well with the number of events required to reject isotropy, 9 
events, in Ref. [10]. 



4 WIMP mass determination 

Several ways of extracting the WIMP mass from direct detection data have been explored. 
The WIMP mass could be extracted from the energy spectrum [24-26] and from the mean 
recoil energy [26] in non-directional detection, from the directional rate as a function of 
9y [11], or from both the angular and energy distributions [27, 28] in directional direct 
detection. The ring angular radius 7 in the differential rate (not integrated over energy) may 
lead to an additional indication of the WIMP mass, since it does not exist for very small 
masses. 

For a given recoil energy and target element, one could obtain some information on the 
mass of the WIMP from the angular radius of the ring. If a ring is present in the differential 
recoil rate at the energy Er, for a target nucleus of mass M one can find m = fiM/(M — ji) 
where 



The m so derived should be the same for all Er (i.e. cos 7 oc \/Er). If a ring does not exist, 

the maximum recoil rate is in the direction of — Vi a b- The ring exists only if the reduced 
mass // is in the range 



Fig. 12 shows a plot of 7 in the differential rate as a function of m (from Eq. 4.1) 
for S, F, and C at three different recoil energies of 5, 20, and 40 keV. The two panels of 
Fig. 12 illustrate the impact of different values of VcaiRot) 312 km/s and 180 km/s, so that 
Mab = 340.2 km/s and 208.8 km/s, respectively. Due to Eq. 2.6, if m > M there is no 
sensitivity of 7 to the WIMP mass m, and 7 reaches a plateau (see Fig. 12). 

It would be difficult to observe a ring with angular radius in the directional rate 7 < 30° 
(the rates at the center and the ring are almost equal), so only for the range of m shown in 
Fig. 13 as a function of Er (between 77130, below which 7 < 30°, and m p i at above which 7 
reaches a plateau), a measurement of 7 could translate into a measurement of m. One cannot 
distinguish m p i at from any mass larger than it, m > m p i at , using the ring radius. Fig. 13 
corresponds to an S detector and VJ a b = 340.2 km/s. Similar ranges can be easily found for 
F and C using Fig. 12. 

Eq. 3.5 can also be used to find the number of events needed to observe the ring in the 
directional differential rate. If we take a narrow enough energy bin of width AEr in keV so 
that the differential rate would be constant in the bin, Eq. 2.12 becomes 




(4.1) 





(4.2) 



K(E,E + AE R ;9 V ) 



dR 



AE R . 



(4.3) 



(IEr d(j)v d cos 6y <t> v 



We choose here a bin AE R ~ 1 keV. For CS 2 , m = 100 GeV/c 2 , E R = 5 keV, Viab = 312 
km/s, and a v = 173 km/s (shown in Fig. 7. a), the exposure needed to observe the ring in the 
differential rate is 4.0 x 10 3 kg-yv-keV / pq^ct^AEr. This corresponds to 19 events needed 
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Figure 12. Ring radius 7, measured in degrees from the direction of — Vi a b; as a function of m for 
three target elements S (black lines), F (green lines), and C (orange lines) and three recoil energies 
E R = 5 kcV (solid lines), 20 kcV (dashed lines), and 40 kcV (dot-dashed lines) for (a) V\ &h = 340.2 
km/s, and (b) Vj a b = 208.8 km/s. Above a certain value of m, 7 reaches a plateau. 
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Figure 13. Sensitivity region where a ring radius can be used to indicated a WIMP mass. The range 
of m is shown between 77130 (black curve) , below which 7 < 30° , and m p i a t (green curve) above which 
7 reaches a plateau (thus m > m p i at ), as a function of Er for VcaiRot = 312 km/s in an S detector. 



in a 1 keV energy bin at 5 keV within 95° of — Vi a b to observe the ring, namely 650 total 
events (over all directions) in the 5-40 keV interval, or 780 events above 20 keV. For the 
same parameters but at = 20 keV, a minimum of 20 events would be required within 64° 
of —Viab to observe the ring which corresponds to an exposure and total number of events 
higher by about a factor of two than those just mentioned. 

We conclude that detecting the ring in the differential rate would require number of 
events at least 6 times larger that those required using the energy integrated rate, i.e. 30 
times more than needed to detect the mean recoil direction according to Ref. [10]. 
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5 Anisotropic logarithmic-ellipsoidal models 



So far we used the isotropic Maxwell-Boltzmann model but here we will show that the ring- 
like feature persists with different characteristics in an anisotropic halo. Various observations 
and numerical simulations suggest that triaxial models with anisotropic velocity distributions 
are a better approximation for galaxy halos [29-31]. In the anisotropic logarithmic-ellipsoidal 
model of Ref. [32], the recoil momentum distribution is anisotropic and can be written as [14] 

$ J± A = i x exp ( [(g/ 2 ^) + q-Viab] 2 \ 

"WV ^{al gq l g + + algq 2 g) X 6XP ^ 2{al g ql g + al gq l g + <Tl g ql g ))> 

(5.1) 

where q xg , q yg , and q zg are components of q in galactic coordinates (see Appendix A), and 
the velocity dispersion matrix in the galactic reference frame is given by, 

'<4g \ 

o* g . (5.2) 
alj 

Since we assume that the WIMPs are on average at rest with respect to the Galaxy, the 
average WIMP velocity — V"i a b in this model is the same as for the 1MB. The equations for 
Qxg-, &ygi and a zg , as well as the constants p and q used to describe the axis ratios of the 
density ellipsoid, and the constant T used to describe the velocity anisotropy, are given in 
Appendix B. 

Eq. 5.1 shows that the angular radius of the ring in this anisotropic model is identical 
to that of the 1MB in Eq. 2.5. The ring is still circular with constant radius at any azimuthal 
angle <py around — Vi a t>, but the value of the rate at the maximum (9y = 7) depends on (fry. 
In the left Mollweide maps of Fig. 14 the ring is almost indistinguishable from the ring in 
the 1MB, but in the right plots of the same figure there is a strong azimuthal variation. 

Keeping p and q constant, the change of the rate at different azimuthal angles is larger 
for larger values of |T|. This can be seen in Fig. 14 where we show a Mollweide map of 
the directional differential recoil rate for m = 100 GeV/c 2 , En = 5 keV in CS2, using the 
Maxwellian distribution in Eq. 5.1 with v esc = 544 km/s and VcaiRot = 312 km/s (on June 2). 
The Solar system is on the major and minor axes in the top and bottom plots, respectively. 
For all four plots, p = 0.9 and q = 0.8. In the top row, T = —1.0 and 16 for the left and 
right plots, respectively, and in the bottom row T = 0.07 and 16 for the left and right plots, 
respectively. The left plots with smaller |T| are similar to those in the 1MB with similar Vj a b 
(Figs. 5.b and 6.c). The right plots with large |T| are very anisotropic with respect to — Vi a b- 
Thus, the ring-like feature can give information on the WIMP velocity distribution. More 
work would be necessary to clarify this issue. 



6 Discussion 



We have pointed out that at low enough recoil energies and for heavy enough WIMPs, the 
maximum of the recoil rate is not in the direction opposite to the velocity of the detector 
with respect to the galaxy, — Vi a b, but on a ring around it at an angle 7 that approaches 90° 
at very low energies. 

For the nuclei of interest in direct detection experiments (at present S, C and F), the 
maximum possible values of 7 are shown in Fig. 3. The ring radius 7 is larger for larger 
values of the average velocity of the WIMPs with respect to the detector. 
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Figure 14. Directional differential recoil rate for m = 100 GeV/c 2 , Er = 5 kcV in CS2, for the 
anisotropic logarithmic-ellipsoidal model with p = 0.9, q = 0.8, v osc = 544 km/s and VcaiRot = 312 
km/s on June 2. Top left (a): Solar system on the major axis and T = —1.0. Top right (b): Solar 
system on the major axis and T = 16. Bottom left (c): Solar system on the minor axis and T = 0.07. 
Bottom right (d):Solar system on the minor axis and T = 16. The color scale/grayscale shown in 
the vertical bars correspond to equal steps between the minimum and maximum values given in units 
of 10~ 6 x (po.3<744/kg-day-keV-sr). The ring is clearly visible and similar to the case of the 1MB in 
(a) and (c). In (b) and (d), |T| is large and the rate is different at different azimuthal angles around 
— Vi a b and is harder to see the ring. 



The ring-like feature can be used as a secondary signature of dark matter since the 
background rate would not have such a feature. The minimum number of events required to 
detect the ring in the energy-integrated rate over relatively low recoil energies is larger than 
required to detect the mean recoil direction, but possible by a factor as low as five, if low 
threshold energies can be achieved in directional detectors. 

In Section 3, for one of the most favorable sets of halo parameters we studied, for a 
100 GeV/c 2 WIMP we found with a simple statistical method that 41 events with energies 
between 5 keV and 40 keV, or 55 events between 20 keV and 40 keV, would be required to 
detect the ring in the energy-integrated directional recoil rate in CS2. With the corresponding 
exposure of 700 kg-yr/po.3C44, there would be 136 events in CS2 above 20 keV (but the ring 
could not be seen if integrating above 40 keV for a 100 GeV/c 2 mass WIMP). 

This number of events necessary to observe the ring, 136, is about 5 times larger than the 
approximately 30 events needed to determine the mean recoil direction according to Ref. [10] 
in the same energy range, with the same WIMP mass and detector material. Using the 
same simple statistical method, we estimated that 13 events would be required to detect the 
forward-backward asymmetry for recoil energies above 20 keV in the same case. This result 
compares well with the approximately 10 events in the same energy range that were found 
in Refs. [9] and [10] to be necessary (using a much more sophisticated statistical method). 

As explained in Section 4, some information can be obtained on the WIMP mass from 
the existence and the angular radius of the ring in the directional differential recoil rate, but 
this would require a number of events at least 6 times larger than those necessary when using 
the energy-integrated directional rate. For example for the most favorable halo parameters 
and a 100 GeV/c 2 WIMP, with the exposure required to observe the ring in a 1 keV energy 
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bin at 5 keV in CS2, 650 events in the 5-40 keV interval, or 780 events above 20 keV would 
be needed. To study the expected precision of the mass determination requires Monte-Carlo 
studies outside the scope of this paper. 

The present threshold of DRIFT is 40 keV. With this threshold, the ring is present in 
S recoils only for WIMPs with m > 600 GeV/c 2 , for which the recoil rates are smaller and 
large number of events would be needed to observe the ring. For a 600 GeV/c 2 WIMP and 
assuming the most favorable halo parameters we studied, 53 events between 40 keV and 50 
keV, which means an exposure of 8.9 x 10 3 kg-yr//9o.3<744, would be required to detect the 
ring in CS2. With this exposure, there would be 315 events above 40 keV. 

We have assumed all along that the recoil directions, including their senses, can be re- 
constructed perfectly in 3 dimensions and there are no backgrounds. Assuming more realistic 
conditions will make the ring-like feature more difficult to observe. 

In conclusion, the ring is the place where most of the lower energy events would come 
from in directional detectors, if the WIMP mass is large enough. This feature can be used 
as a secondary signature for dark matter, given that it requires more events to be detected 
than other signatures already discussed in the literature. The detailed features of the ring 
depend on the WIMP mass and velocity distribution. Further study is needed to find the 
optimal strategy and the best achievable precision in the determination of the WIMP mass 
using the ring. The ring-like feature could also be used to test the shape of the WIMP 
velocity distribution. For example, we showed that logarithmic-ellipsoidal distributions give 
an azimuthal pattern around the ring, while the ring remains circular. Further study is 
required to clarify this issue too. 
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A Complete transformation equations from the detector frame to the 
Galactic reference frame 

The presence of q • V] a b in /i a b means that in order to compute the differential rate we 
need to orient the nuclear recoil direction q with respect to Vi a b- We present here complete 
transformation equations for q and Vi a b to go from the detector frame to the Galactic frame. 

We define a reference frame fixed to the laboratory and orient its axes so that the xy 
plane is horizontal, the x-axis points North, the y-axis points West, and the z-axis points to 
the Zenith. We denote its unit coordinate vectors as Jv, W and Z, respectively. The detector 
is at some orientation in the laboratory. We define the detector frame with X, Y, Z cartesian 
axes fixed to the detector. The unit coordinate vectors of the detector frame are X, Y and 
Z. The transformation between the lab frame and the detector frame is given in Eqs. Al to 
A3 of Ref. [13] in terms of "direction cosines" . q is given in the detector reference frame as 
q = qx X + qy Y + qz Z, and we can write it in the lab frame using Eq. Al of Ref. [13], 
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q w 



q = q n N + q w VV + q z Z, where 

qxax + qyOLY + gz«z, 
: qxfix + qvPv + qzfiz, 
q z = qxix + qviv + qziz, (A.l) 

and cti, Pi and 7$ are the direction cosines between the two sets of cartesian coordinates of the 
lab and detector frames, for i = X, Y, Z. For example ay is the cosine of the angle between 
the N and Y directions, and fiz is the cosine of the angle between the W and Z directions. 

We would like to write q in the Galactic reference frame. Using Eq. A8 of Ref. [13], q 
can be written in the equatorial frame, 



q = [-q n sin(Ai a b) cos^b) + q w sin(i° ab ) + q z cos(Ai ab ) cos(i° ab )] x e 
+ [-q n sin(Ai a b) sin(i lab ) - q w cos(i° ab ) + q z cos(Ai ab ) sin(t° ab )] y e 
+ [q n cos(Ai a b) + q z sin(Ai ab )] z e , 



(A.2) 



where x e , y e , and z e are the unit coordinate vectors of the geocentric equatorial inertial 
(GEI) frame: its origin is at the center of the Earth, its x e -axis points in the direction of the 
vernal equinox, its y e -axis points to the point on the celestial equator with right ascension 
90° (so that the cartesian frame is right-handed), and its z e -axis points to the north celestial 
pole. t° ab is the laboratory Local Apparent Sidereal Time in degrees, and Ai ab is the latitude 
of the lab. 

The transformation from the Galactic frame to the equatorial frame is given by 



(A.3) 
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where 




(A.4) 



and x g , y g , and z g are the unit vectors of the Galactic reference frame. Recall the definition 
of the Galactic coordinate system: its origin is at the position of the Sun, its x g -axis points 
towards the Galactic center, its y s -axis points in the direction of the Galactic rotation, 
and its z g -axis points to the north Galactic pole. Note that these coordinates are related 
to Galactic longitude I and latitude b by (x g , y g , z g )=(cos 6cos I, cos 6sin I, sin b). We have 
q = q xg x 9 + q yg y g + q zg z g , where 



qxg — qn 
+ q z 

Qyg = In 
+ c 

q Z g — qn 
+ q z 



- [a x cos(t lab ) + a y sin(t lab )] sin(Ai ab ) + a z cos(Ai ab )) + q w (a x sin(t lab ) - a y cos(t lal 
[a x cos(t! ab ) + a y sin(t lab )] cos(Ai ab ) + a z sin(Ai ab )) , 

■ [b x cos(t lab ) + by sin(£ lab )] sin(Ai ab ) + b z cos(Ai ab )^ + q w (b x sin(£ lab ) - b y cos(t° ab ' 
[b x cos(t° ab ) + by sin(^° ab )] cos(Ai ab ) + b z sin(Ai ab )^) , 

- [c x cos(t° ab ) + Cy sin(t° ab )] sin(Ai ab ) + c z cos(Ai ab )^ + q w (c x sin(t° ab ) - c y cos(t lab ; 
[c x cos(t lab ) + Cy sin(t lab )] cos(Ai ab ) + c z sin(Ai ab )J . (A.5) 
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As done in Ref. [33], the transformation matrix Aq can be found using the definition 
of the Galactic pole and center in the International Celestial Reference System (ICRS). The 
north Galactic pole can be defined by right ascension ac = 192°. 85948 and declination 
5q = +27°. 12825 in the ICRS. The origin of Galactic longitude is defined by the Galactic 
longitude of the ascending node of the Galactic plane on the equator of ICRS, which is 
la = 32°. 93192. The angles a<j, Sq and la should be regarded as exact quantities, and they 
can be used to compute the transformation matrix Aq [33], 

/ -0.0548755604 +0.4941094279 -0.8676661490 \ 
A G = -0.8734370902 -0.4448296300 -0.1980763734 . (A.6) 
\ -0.4838350155 +0.7469822445 +0.4559837762 / 

The direction of q in the Galactic rest frame can be specified by the Galactic longitude 
(lq) and latitude (b q ) of q, 

q xg = cos b q cos l q , q yg = cos b q sin l q , q zg = sin b q , (A.7) 

where q xg , q yg , and q zg are given in terms of q n , q w , and q z in Eq. A. 5. 

The velocity of the lab with respect to the center of the Galaxy can be divided into 
four components: the Galactic rotation velocity VcaiRot at the position of the Sun (or Local 
Standard of Rest (LSR) velocity), Sun's peculiar velocity Vsoiar in the LSR, Earth's transla- 
tional velocity VgarthRev with respect to the Sun, and the velocity of Earth's rotation around 
itself VEarthRot- We also need to write each component of Vi a b in the Galactic frame and 
compute, 

q • V lab = q • VcalRot + q • V So ia r + q • V Ea rthRev + q ' V Eart hRot- (A. 8) 

A.l Galactic rotation 

The velocity of the Galactic rotation VoaiRot is defined in the Galactic reference frame, 

VcalRot = ^GalRotJV (A. 9) 

Recall that VcaiRot is the Galactic rotation speed (i.e. the local circular speed), and y g is in 
the direction of the Galactic rotation. Thus we have 

q • V Ga lRot = gj/gVbalRot (A. 10) 

where q yg is given in Eq. A. 5. 

The standard value of VcaiRot for the 1MB is 220 km/s [34]. As discussed in Ref. [10], 
recent studies have found other values for VbalRot • O ne analysis which was argued that it has 
used overly restrictive models found VcaiRot = (254 + 16) km/s [35]. Another study found 
^GaiRot — (236 ± 11) km/s assuming a flat rotation curve [36], while a different study found 
values ranging from VoalRot = (200 ± 20) km/s to VcaiRot = (279 ± 33) km/s depending on 
the model used for the rotation curve [37]. We take VcaiRot = 180 km/s and 312 km/s as 
low and high estimates. 

A. 2 Solar motion 

The Sun's peculiar velocity in the LSR is, 

Vsoiar = Ui g + Vy g + Wi g . (A. 11) 
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Thus, we can compute q • V"s iar as 



q • V Solar = q xg U + q yg V + q zg W. (A. 12) 

Ref. [38] has re-evaluated the velocity of the Sun with respect to the LSR and found 
that the classical determination of U, V, and W is undermined by the metallicity gradient 
in the disc. In their study they find (U,V,W) & = (ll.l±g;f|, 12.24±g;g, 7.25±° f 6 ) km/s, 
with additional systematic uncertainties ~ (1,2,0.5) km/s. The new values are extremely 
insensitive to the metallicity gradient within the disc. 

A. 3 Earth's revolution 

The velocity of the Earth's revolution around the sun is given in terms of the Sun ecliptic 
longitude X(t) as [39] 

V E arthRcv = V©(A) cos /3 X sin (A - A x )x g + cos (3 y sin(A - X y )y g + cos (3 Z sin(A - X z )z g 

(A.13) 

where V® = 29.8 km/s is the orbital speed of the Earth, V©(A) = V©[1 — esin(A — Ao)], 
e = 0.016722, and Ao = 13° ± 1° are the ellipticity of the Earth's orbit and the eclip- 
tic longitude of the orbit's minor axis, respectively, and = (—5°. 5303, 59°. 575, 29°. 812) 
and Xi = (266°. 141, -13° .3485, 179°. 3212) are the ecliptic latitudes and longitudes of the 
(x g ,y g ,z g ) axes, respectively. 

The Sun's ecliptic longitude A can be expressed as (Ref. [40] and p. 77 of Ref. [41]), 

A = L + (1°.915 - 0°.0048T ) sin g + 0°.020 sin 2g, (A.14) 

where L = 281°. 0298 + 36000°. 77T + 0°.04107UT is the mean longitude of the Sun corrected 
for aberration, g = 357°. 9258 + 35999°. 05T + 0°.04107UT is the mean anomaly (polar angle 
of orbit). UT is the Universal Time in hours, and To is the time in Julian centuries (36525 
days) from 12:00 UT on 1 January 2010 to the previous midnight. 
Thus, we can compute q • VgarthRev as 

q • V EarthRcv = U®(A) cos(3 x sin(A - X x )q xg + cos/3 J/ sin(A - X y )q yg + cos f3 z sin(A - X z )q zg 

(A.1'5) 

A. 4 Earth's rotation 

Finally, we want to compute VgarthRot! the velocity of Earth's rotation around itself. We 
have 

V E arthRot = -U Rot Eq COS Ai ab VV, (A.16) 

where VkotEq is the Earth's rotation speed at the equator, and is defined as Ua ot Eq = 
27T.R®/(1 sidereal day). The Earth's equatorial radius is R@ = 6378.137 km, and one sidereal 
day is 23.9344696 hr= 86164 s. therefore VaotEq = 0.465102 km/s. 

We use Eq. A8 of Ref. [13] to write VV in terms of the equatorial frame coordinates, 

V Ea rthRot = -^RotEq cos Ai ab ( sin(t lab )x e - cos(t° ab ) y^j , (A.17) 
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Using the the Galactic to equatorial transformation, we have 

V E arthRot = -VR ot Eq COS A lab | (a x Sm(t^ h ) - d y COs(t° ab )J X g 

+ (b x sin(t 1 ° ab ) - b y cos(t^ h )jy g + (cx sm(t^ h ) - Cycos^))^}. (A.18) 
Thus we can compute q • VEarthRot as 

q • V EarthRot = -^RotEq cos Ai ab | (a x sin(t° ab ) - a y cos(t° ah )jq xg 

+ (b x sin(i° ab ) - by cos(i lab )J q yg + (c x sin(i° ab ) - c y cos(i lab )J q^. (A. 19) 
Inserting q xg , q yg and q zg from Eq. A. 5, we find 



q • V Eart hRot = -VRotEq cos Ai ab <! ( [at + b y + c? - a x - b x - c x ] cos(tJ ab ) sin(i° ab ) sin A L 



a.b 



+ [a x a y + b x b y + c x c y ] (cos 2 (^ ab ) - sin 2 (i lab )) sinAi ab 

+ [a x a z + b x b z + c x c z ] sin(i° ab ) cos Ai ab + [a y a z + b y b z + c y c z ] cos(i° ab ) cos Ai ab ) q v 
+ ([«x + b 2 x + c 2 ] sin 2 (t 1 ° ab ) + [a 2 + b\ + c 2 ] cos 2 (^ ab ) 
- 2[a x a y + + c x c j/] cos(i[ ab ) sin(t[ ab ) )q w 



+ [K + b x + c x - a - b - c] cos(t° ab ) sin(t 1 ° ab ) cos A 



lab 



lab 



+ [a x a y + b x b y + c x c y ] (sin 2 (t lab ) - cos 2 (i lab )) cos A 

+ [a x a z + b x b z + c x c z ] sin(i lab ) sin Ai ab - [a y a z + b y b z + c y c z ] cos(^ ab ) sin Ai ab ^ 2 |. 

(A.20) 

From Eq. A. 6, we have aiCLj + bibj + CjCj = 5ij where i,j = x,y,z. Therefore, Eq. A.20 will 
be simplified to 

q • V EarthRot = -q w VRotEq cos Ai ab . (A.21) 

A. 5 Total Velocity 

Using Eqs. A. 10, A. 12, A. 15 and A.21 we can compute q • Vi ab as 



q • Vi ab = q xg \U + V©(A) cos /3 X sin(A - X x )j + q yg \V GalRot + V + V® (A) cos /3 y sin(A - A,,) 
+ 1zg(w + V©(A) cos fi z sin(A - A z )) - q w V Rot Eq cos A [ab , (A. 22) 

where q xg , q yg , and q zg can be written in terms of q n , q w , and q z using Eq. A. 5. We can also 
write q • Vi ab in terms of only q xg , q yg , and q zg using Eqs. A. 10, A. 12, A. 15 and A. 19, 



q • Vi ab = q xg \U + Ve(A) cos /3 X sin(A - X x ) - VRotEq cos A iab (a x sin(i° ab ) - a y cos(t lab )) 

+ Qyg ( VbaiRot + V + V©(A) cos f3 y sin(A - X y ) - V Rot E q cos Ai ab (b x sin(t lab ) - b y cos(^ ab )) 



+ qzg[ W + V©(A) cos p z sin(A - \ z ) - VRotEq cos Ai ab (c x sin(f lab ) - c y cos(i lab )) j . 

(A.23) 
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B Logarithmic-ellipsoidal models 



In one logarithmic-ellipsoidal model provided by Evans, Carollo, and de Zeeuw [32], the Solar 
System is on the long (major) axis of the halo density ellipsoid, and 

V 2 

2 ^GalRot (-n -1 \ 

X9 (2 + r)(p- 2 + q- 2 -l)' [ 1 



_2 _ ^GalRotV"^ 21 f-o n\ 

°va ~ nt 2 , 2 TV' \ D - Z ) 



VGalRot(2p- 2 ~ 1 
" 2(p" 2 + q-2-l) ' 

yGalRot(2q" 2 -l 

Z9 2(p- 2 + q- 2 - 1) 



'2 _ • GalRotV"^ -V (x> q\ 

a zo ~ o^-2 i —2 TV' 



In another model by the same authors, the Solar System is on the intermediate (minor) axis 
of the halo density ellipsoid, and 



2 ^GalRotP 4 

X9 (2 + T)(l + q-2-p-2V 



^alRot(2-p- 2 ) 

2(l + q- 2 -p- 2 , 

^aiRot(2q- 2 -P^ 
23 2(l + q-2-p-2) 



_2 _ r GalRotV"l r I fj3 f?\ 



p and q are constants used to describe the axis ratios of the density ellipsoid, and T is a 
constant used to describe the velocity anisotropy. 

We have studied 12 different combinations of parameters in the logarithmic-ellipsoidal 
model: 6 combinations for the Solar system on the major axis of the halo density ellipsoid 
and 6 on the minor axis. Eight of the cases are taken from Table II of Ref. [14] and four from 
Table II of Ref. [42]. The combinations of parameters for the Solar system on the major axis 
are 



p 


= 0.9, 


q = 


0.8, 


r = 


-1.78, 


P 


= 0.9, 


q = 


0.8, 


r = 


16, 


p 


= 0.9, 


q = 


0.8, 


r = 


-1.0, 


p 


= 0.9, 


q = 


0.8, 


r = 


-1.33, 


p 


= 0.72. 


i q = 


= 0.7. 


, r = 


= -1.39 


p 


= 0.72. 


- q = 


= 0.7. 


, r = 


= -1.6, 



(B.7) 

and the combinations of parameters for the Solar system on the minor axis are 

p = 0.9, q = 0.8, T = -1.78, 

p = 0.9, q = 0.8, r = 16, 

p = 0.9, q = 0.8, T = 0.07, 

p = 0.9, q = 0.8, r = -0.62, 

p = 0.72, q = 0.7, r = 4.02, 

p = 0.72, q = 0.7, T = 2.01. (B.8) 
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